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ERROR  THEOKf  OF  INTERSECTION  PHOTOGRAMMETRI 


ABSTRACT 


Formulas  are  developed  for  evaluating  the  propagation  of  the 
various  errors  present  in  a  system  of  photogrammetric  cameras  in 
triangulating  the  space  position  of  a  photographed  point.. 
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I,  INTRODUCTION 


If  at  different  points  oil  the  surface  of  the  earth  photagrammetric 
cameras  are  established  in.  order  to  record  points,  and  if  the  orientation 
of  the  photographs  is  known  and  the  corresponding  plate  coordinates  are 
measured,  then  the  spatial  position  of  such  points  may  be  computed. 

Besides  the  errors  in  the  geodetic  reference  data  of  the  stations,  the 
errors  of  the  parameters  of  orientation,  and  of  the  measured  plate  coordi¬ 
nates  contribute  to  errors  in'  the  computation  of  the  individual  spatial 
direction  to  any  target  point,  and  consequently  lead  to  errors  of  the 
spatial  coordinates  of  the  targeto  This  paper  deals  with  the  determi¬ 
nation  of  the  influence  of  the  individual  residuals  and  gives  formulas 
for  the  combined  effect  of  these  errors. 

A  study  of  the  error  theory  of  intersection  photogrammetry  must 
follow  steps  which  are  typical  of  the  application  of  the  method.  A 
basic  characteristic  of  the  triangulation  method  is  the  fact  that  each 
camera  in  the  course  of '"the  reductions  is -treated  as  an  independent 
unit  up  to  the  moment  when  the  direction  to  a  target  is  obtained.  Such 
a  direction  may  be  expressed  by  two  position  angles  (e.g.,  Azimuth  and 
Zenith  distance).  Only  in  the  final  step  in  the  triangulation  procedure 
will  the  directions  from  the  individual  cameras  be  combined  for  the  deter¬ 
mination  of  the  most  probable  point  of  intersection.  Consequently,  the- 
problem  is  to  develop  first  the  error  theory  of  an  individual  photogram- 
. metric  camera,  and?  later  to  determine  the  propagation  of  the  individual 
station  errors  in  the  triangulation  procedure. 

II.  THE  .ERROR  THE OKI  OF  A  PHOTOGRAMMETRIC  CAMERA 

1.  The  orientation  of  a  photogrammetric  camera 

•  a- 

We  introduce  a  rectangular  cartesian  coordinate  system  in  such  a 
way  that  its  origin  is  the  center  of  projection  and  its  orientation  is 
consistent  with  the  local  geodetfc  reference  system.  Thus,  the  X,  Y, 
plane  represents  the  horizontal  plane  at  the  station.  .  The  X-axis  may 
be  chosen  positive  to  t he  north.  The  Z-axis  is  perpendicular  to  the 
X,  Y  axes,  positive  to  the  zenith.  Further,  we  always  consider  the 
interior  and"  exterior  orientation  of  a  plate  simultaneously;  that  is, 
we  consider  the  plate  in  its  spatial  position.  Consequently,  the  plate 
coordinate  system  -  -primarily  a  plane  system  -  will  be  introduced  as 
another  rectangular  cartesian  spatial  system  with  its  origin  again  at 
the  center  of  projection.  The  plate  is  introduced.as  a  diapositive. 

Figure  1  shows  the  situation. 
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L.  is  the  center  of  projection .  Its  geodetic  coordinates  denoted 
by  xL,  yL  and  zL  may  be  expressed  in  any  suitable  rectangular  coordi¬ 
nate.  system  which  need  not  necessarily  be  a  cartesian  system.  In  L. 
we  establish  two  rectangular  cartesian  systems.  The  X-Y-Z  system, 
denoted  as  the  station  system,  is  expressed  by  the  unit  vectors  i, 
j,  K.  In  this  system,  a  spatial  target  point  R  is  determined  by  the 
coordinates  X,  Y,  Z.  Its  image  point  is  r,  denoted  by  the  coordi¬ 
nates  u,  v,  w  respectively.  The  system  is  oriented  .according  to  the 
geodetic  reference  system.  The  other  cartesian  system  expressed  by 
the  unit  vectors  x,  j,  fc  denotes  the  camera  system.  Its  orientation 
relative  to  the  station  system  is  expressed  by  the  three  elements  of 
exterior  orientation  denoted  by  the  azimuth  angle  A  (counted  clock¬ 
wise  from  the  south),  the  tilt  angle  v  of  the  plate  perpendicular, 
and  by  the  swing  angle  /f  of  the  fiducial  mark  system.  The  length  of 
the  plate  perpendicular  -  the  principal.,  distance  -  is  denoted  by  d. 

A  x  and  Ay  are  the  coordinates  of  the  .principal  point  P  in  any 
rectangular  plane  plate  coordinate  system  established  by  fiducial 
marks  and  denoted  by  x  and  y„  x  and  y  are  the  plane  plate  coordinates 
of  the  image  point  r  in  the^briented  x,  y  system. 

Further,  we  introduce  the  image  vector  7  and  the  object  vector  7 
respectively. 

Then? 

7  =■  iu  +  jv  +  kw  =  +id  -  $c  ”  (1) 

and  7  =■  iX  +  jY  +  kZ  . 

The  transformation  of  the  triple  vector  d),  3,  ^  into  the  triple  vector 
i,  j,  k  is  obtained  from  the  transformation  matrix 


fhe  individual  transformation  matrices  are: 

1. )  for  a  rotation  in  aiimuth  denoted  temporarily  by#: 
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with  fortnulas  (1)  and  (2),  we  obtain? 

A  * 

+id  =  -id  cos  A-  sin  v  -  jd  sin  A  sin  v  +  kd  cos  v 
-jx  =  -ix  sin  A  +  jx  cos  A 

-ky  *  -iy  cos  A  cos  v  -  jy  sin._A  cos  v  -  ky  sin  v 


r  =  i(-x  sin  A  -  y  cos  A  cos  v  -  d  cos  A  sin  v)  +  j(+x  cos  A  -  y  sin  A  cos  v  — 

d  sin  A  sin  v)  +  k(-y  sin  v  +  d  cos  v) 

and  again  from  formula  (1) 

u  =  -x  sin  A  -  y  cos  A  cos  v  -  d  cos  A  sin  v 

v  ■  +x  cos  A  -  y  sin  A  cos  v  -  d  sin  A  sin  v  (3) 

w  =  -y  sin  v  +  d  cos  v 


— * 

The  coordinates  X  and  Y.  of  a  point  along. the  vector  R  for  a  certain 
elevation  Z  are' therefore: 


X  =  -  (+x  sin  A  *  y  cos  A  cos  v  +  d  cos  A  sin  v) 
+y  sin  v  -  d  cos  v 

Y  =  Z  (~x  co3  A  +  y  3in  A  cos  v  +  d  sin  A  sin  v) 
+y  sin  v  -  d  nos  v 


From  Figure  1  we  obtains 

_  y  t  |  f 

x  =  (x  -  ax)  cos  <  ~  (y  -&y)  sin  tc 

(5) 

y  e  (y  -  Ay)  cos  k  +  (x  -  Ax)  sin  *c 


Substituting  these  values,  in  (ii),  we  have: 

Z  (|[(y-Ay)cos  K+  (x“Ax)  sin.tfj  cos  v+d  sin  vj  cos  A+  £(x-ax)cos<  -(y-Ay)  sin^cjsin 

O  -ay)  cosK  +  (x-dx)  sin  kJ  sin  v-d  cos  v  (<: 


*These  formulas,  are  in  agreement  with  formulas  Mo.  12  derived  by 
0.  r.  Gruber  .in  Ferienkufs  in  Photogrammetrle,  Page  27 • 
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«4 


y  =  2  ^[(y~Ay)C03'<  +  (^~-Ax)s3jiaQC0S  v+d  sin  vjsifa  A-f(x~£x_)  cosx  -(y-Ay)sin/fjcos  cont^ 


£(y-Ay)  cos*  +(x”Ax)sln>cJ-siii  v-d  cos  v 


») 


With  Z  =  1,  these  formulas  reduce  to  formulas  10a  for  the  standard 
coordinates  £  ,  ij  obtained  in  BRL  Report  No.  7'8U  ■* 

Sometimes-,  especially  in  aerial. photogramraetry,  it  becomes  necessary 
to  express  the  direction  of  the  camera  axis,  by  two  tilt  angles.  Let  us 
denote  the  tilt  in  the  x-direction  by  a  and  the  tilt  in  the  y-direction 
by  to,  as  in  Figure  2.  We  have: 

~t  =  iu  +  jv  +  kw  =  id  +  5y  -  fex  .r  •.  (7) 


The  individual  transformation  matrices  are: 
l.}for  a  rotation  in  a-tilt,  denoted  temporarily  by* 


i*  i*  -  i  cos  co  -  -5*  sin  co 


*H.  Schmid  -  ‘'Spatial  Triangulat'xon  by  Means,  of  Phot  Pgr  amme  try." 
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and  the  combined  transformation  matrix  is: 


i  »  i  sin.  a  cos  co  +  j  sin  co  +  k  cos  a  cos  co; 

j*=-i  sin  a  sin  co  +  j  cos  co  -•  k  sin  co  cos  a 

£  =-i  cos  a  +  k  sin  a  (8) 

with  formulas  (7)  and  (8)  we  obtain: 

+id  =  +id  sin  a  cos  co  +  jd  sin  co.  +  kd  cos  a  cos  co 

3  -iy-  sin  a  sin  co  +  jy  cos  co  -  ky  cos  a  sin  ai 

4ac  =  cos  a  -lac  siii  a  x 

7  =  i(d  sin  a  cos  co  -  y  sin  a  sin  co  +  x  cos  a) 

+j(d  sin  co  +  y  cos  co)  +  k(d  cos  a  cos  co  -  y  cos  a  sin  co  -  x  sin  a) 

i 

and, again  with  formula  (7) 

u  =  d  sin  a  cos  co  -  y  sin  a  sin  co  +  x  cos  a 
v  3  d  sin  co  +  y  cos  co  - 

w  =  d 'cos  a  cos  co  -  y  cos  a  sin  co  -  x  sin  a  ( 9 ) 

The  coordinates  X  and  Y  for  a  point  along  the  vector  R  for  a  certain 
elevation  2  are  therefore: 
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(10) 


X  _  Z(d  sin  g  cos  co  -  y  sin  a  sin  u>  +  x  CQ3  a) 
d  cos  a  cos  a>  -  y  cos  a  sin  co  -  x  sin  a  , 

j  =  Z.  (d  sin  m  +  y  cos  to)  ■  _ _ _ 

d  cos  a  cos  go  -  y  cos  a  sin  a>  -  x  sin  a 

From  Figure  2,  we  obtain  y  ==  (y-Ay)  cos  x  -  (x-£x)  sin'/c 

x  =  (J-iSy)  sinx  +  (x-Ax)  cds  K 

Substituting  formula  (11)  into  formula  (10) 


(11) 


-  zff(x-Ax)  cosCKy.-&y)sinjflcos  a+‘;d  cos .  co-E(y-Ay)  cos*-(x-Ag)sin^lsin  coTsin  a) 
~  £d  cos  o>"'[r(y”Ay)cos#C“(x”Ax)sinx^in  oojcos  a-^x-Ax)  cosK+  (y-Ay )  sin >c]din  a . 

(12)  * 

y  b  dll  y-Ay)  cos>C-(x-Ax)  sin^jcos  go  +  d  sin  go}-  _ 

{d  cosa>^(y-Ay) cos K~ (x-Apc) djj)^m»JcQSC^x-Ax) cos^+  (y-^y )sin xjsin  a' 

For  Z  =  1  (unit  distance),  formulas  (12)  reduce  to  the  expressions  given 
in  the  BRL  Report  Wo.  78iv*fc  for  the  standard  coordinates  ^  and  >?  . 

2 «  Differential  Formulas 

We  will,  now  consider  the  differential  changes  on  the  image  vector 

"r*  and  object  vector  R  respectively.  These  changes,  are  caused  by  dif¬ 
ferential  changes  of  the  individual  parameters  of  the  solution. 

We  denote  the  spatial  position  of  the  optical  axis  with  respect  to 
the  station  system  either  by  the  azimuth  A  and  the  zenith  distance  v,  or 
by  the  a-tilt  and  oo-tilt  angles.  The  .corresponding  differential  changes 
are  dA,  dv,  and  da,  da  respectively.  The  swing  angle  of  the  fiducial 
mark  system,  is  denoted  by  tC  and  its  differential,  change  by  d  fC  .  The 
principal  distance  of  the  camera  is  denoted  by  d  and  the  coordinates  of 
the  origin  of  the  fiducial  mark  system  denoted  by’  0  are,  in  the  oriented 
x3y~system,  xq,  and  yQ«  The  corresponding  differential,  changes  are  dd, 

dx  and.dy  respectively^  The, location  of  any  target  point  on  the  plate  is 
determined  by  a. fiducial,  mark  system  and  expressed  by  x  and  yj  the  dif¬ 
ferential  changes  being  dx  and  dy.  Finally,  the  spatial  location  of  the 
center  of  projection  (L)  within  a  certain  geodetic  reference  system  is 

^Formulas  12  are  in  agreement  with  formulas  13  derived  by  0.  v.  Gruber. 
Compare  footnote  on  .page  *7  « 

#*cf  reference  on  page  3  ° 

HiurDIJT".! 

•■vrvyy-; 
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ELEMENTS  OF  ORIENTATION 


given  by  the  coordinates  xL,  and  zu  and  the  corresponding  differential 
changes  arc  dxL,  dyt  and  dz^  » 

Coordinates  in  the  various  systems  according  to  this  nomenclature  are 
summarized  in  the  following  tables 


Azimuth,  zenith 
distance-system 


a-tilt,  <o-tilt- 
_ system _ 


tv 


O  O 
CO  -P  ‘H 
O  G  -P 
-P  0)  o 
g  O 

•H  fl>  O 
T3  G  ^ 

c.  -p  a 

o 

0,<H  Vh 

o  o  o 


<u 

5 


*h  «J. 

•h  n  nj 
■p  a>  o 

*H  (— I  "H 
CO  WM->' 
O  G  O. 

Pu  ro  o  i 


x  -  coordinate 
y  -  coordinate 
z  -  coordinate 


xL  -  dx^ 

yL  1  *y-L 

z.  -  dz, 


Azimuth 

Zenith  distance 
Swing 


A  -  dA 


v  -  dv 


#  *  dK 


xL  -  dxL 

5! i.  1  dyL 

zL  -  dzL 


a-tilt  a.  -  da. 
co-tilt  co  -  do 
K  -  dAC 


‘  *=• 

0  O  M 

.  u 

r  p 

»  33  « 

co  ■  e 

o 

■H 

u 

g'Syn 

G  .  ct)  _ 

G 

•rl  G  »H  E- 
■p  tH  o  O’ 

C.  bn. 3  -p 

M 

O  -H  t)  p 

O.  Ch  -f-1  Pa 

O  O  ^  (0 

Principal.,  .distance 
X  -  coordinate 
y  -  coordinate 


d.  -  dd 


xQ  i  dx 

1  dy 


d  *  dd 


x  -  dx 
o 

y^  -  dy 


63 


so 


3  s 

Pu  S 


8 

CJ> 


1,5  - 

>H  £ 

O  <L' 

3  -P 

t)  tj  n 

-H 

(*<  E  co 


x  -  coordinate 
y  -  coordinate 


x  -  dx 

—  +  ~ 
y  ~  dy 


x  -  dx 
y  -  dy 


11s 


A  differential.,  translation  ds  (dzj_ ,  dyL  and  dz^ )  of  the  center  of 

projection  causes  a  change  in  the  vector  r,  corresponding  . to  the 
relation  ds  =  idx^_  +  jdy^  +  kdz^  ( 

The  three  differential  rotations  dA?  dv  and  d  K  cause  a  differen¬ 
tial  rotation  da  on  the  vector  r  .  From  Figure  1  and  31  we  obtain 

*  da  =  +i  (-d  sin  v  cos  A  +  dv  sin  A) 

+j  (-d>C  sin  v  sin  A  -  dv  cos  A) 

■*7  +k  (+d  K  cos  v  +  dA)  .  ^ 


dZs/nj  cos  A 

dK  s/n  V  sin  A 


or/  <K 


4 


\v  I 

\  I 

Y  I 
V 


Figure  3 


and  the  corresponding  differential  rotation  in  the  tilt-systera  is 
with  Fig®  2  and 


da  =  +i(-d  tc  cos  01  sin  a  -  d  01  cos  a) 
+j(~d  k  sin  01  +  d.  a) 

+k(-d  tc  cos  oi  cos  a  +  d  a>  sin  a) 


(  "r"  +  )  =>  i(u  +  §xi )  +  j(y  +  £y)  +  k(w  +  fi'W) 

\ 

=  i(d  +  dd)  -  ^(x  +  dx)  ~  1c(y  +  dy) 


(16) 


(15). 
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and  from  (3) 

(  r  +  dr  )  -  i(u-dd  cos  A  sin.  v  -  dx. sin  A  -  dy  cos  A  cos  v) 
+j(v-dd  sin  A  sin  v  +  dx  cos  A  -  dy  sin.  A  cos  v) 
+k(w+dd  cos  v  -  dy  sin  v) 


(17) 


and  the  corresponding  computations  with  formulas.  7  and  9  give  for  the 
tilt-system: 

(  *r  +  dr^ )  =  i(u  +  £u)  +j(v  +  Sv)  +  k(w  +  8w)  =  i(d+dd)+j‘(y+dy)-lc(x+dx) 

(18)  • 

(  r  +  dr  )  =  i(>»+dd  sin  a  cos  oo  +  dx.  cos  a  -  dy  sin  a  sin  <o) 

+j(v+dd  sin  co  +  dy  cos  to)  (19) 

+k(w+dd  cos  a  cos  co  -  dx  sin.  a  -  dy  cos  a  sin  co) 


The  differential  changes  dx  and  d(y  are  two  further  translations  in 
the  i,  j,  1c  system.  With  formulas  (5)  and  (11)  respectively,  we  have: 

dx  “  dy  sin  K  +  dx  cos  k 

dy  =  dy  cos  *c  -  dx  sin  k 

and 


dx  =  dx  coa  K  -  dy  sin  fc 
dy  =  dx  sin  K  +  dy  cos  /c 


and  therefore,  with  (17)  and  (19)  respectively 
d?  =  i[^dx(+sin  A  cos cos  A  cos  v  sin/f)+dy(sin  A  sin/C-cos  A  cos  v  cos  *c)^ 
+j  l+dx(+cos  A  cos sin  A  cos.  v  sin < )-dy(+cos  A  sin<+sin  A  cos  v  cos<)  [ 

r  -  -  -l  (20) 

+k/">d?c  sin  v  sin  K  -  dy*sin  v  cos  tc  I 
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and 

d? *»  i[dx(cos  a  cos/C+sin  a  sin  co  sin  >c)+dy(cos  a  sin/c~sin  a  sin  co  cos /cfj 
j|-dx  cos  co  sirr<+dy  cos  co  cos/c/ 

r-  -l(2l) 

Wdx(“3in  a  cos/C+cos  a  sin  co  sin /C)“dy(sin  a  sinx+cos  a  sin  co  cos  «)1 

The  resulting  differential  change  of  the  vector  r  is  obtained  by 
superimposing  the  individual  changes,,  .For  A,  v-system,  we  have  with 
formulas . (13),  (lii),  (17),  and  (20). § 

\ 

dr  =  £da'i^J*-i  [~dd  cos  A  sin  v~dx  sin  A-dy  cos  A  cos  v 

“dx(sin'  A  cos/C+cos  A  cos  v  sin  x)+dy(sin  A  sinX-cos  A  cos  v  cos /c) 

+dxJ 

.+ j  J-dd  sin  A  sin  v+dx  cos  A-dy  sin  A  cos  v+dx(cos  A  cos  tC -sin  A  cos  v  sin/^ 

-dy(cos  A  sin/c+sin  A  cos  v  cos  <“)+dyL^  (22) 

+k  fa*  cos  v-dy  sin  v-dx  sin  v  sinX“dy  sin  v  cos<C  +  dz^J 

and  for  the  a-co-tilt  system  with  formulas  (13),  (1E>),  (19)  and  (21), 
dr  =  t  da"^j+i^-dd  sin  a  cos  co+dx  cos  a-dy  sin  a  sin  co+dx(cos  a  cos  x+  sin  a  siiiXsinco) 
+dy(cos  a  sin  X -sin  a  cos/csin  co)+dXj_J 
+j[dd  sin  co+dy  cos  w-dx  cos  co  sinX+dy  cos  co  cos/c+  dyLJ 

+kjdd  cos  a  cos  co-dx  sin.  a-dy  cos  a  sin  co+dx(-3in  a  cosX+cos  a  sin  co  sin*) 

.  “dy(sin  a  sinX.+  cos  a  sin  co  cos  xO+dz^J  (23) 

,  1 
The  vector  product  ^do'^r’J in  formulas  (22)  and  (23)  may  be  written 
with  formulas  (Hi)  and“(l5)  respectively,  ass 


Tdo  r*l  = 

i 

-d  /f  sin  v  cos  A+dv  sin  A 

j  ■ 

-d  X  sin  v  sin  A-dv  cos  A 

k  1 

+d/Ccos  v+dA  1 

L  J 

u 

V 

w  | 
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and  in  the  tilt  systems 

r  1  j  k 

Idcr  r  1=  -dtfcos  oi  ,sin  a-doj  cos  a  -d)£sin  co+da  -dtfcos  05  cos  a+du  sin  a 

u  T  M  '  (25) 

A  reduction  of  -these  determinants  gives  for  (2l±) 
jdcrr*]=  i[jd/c(w  sin  v  sin  A+v  cos  v)-dv  •  w  cos  A-dA  *  v] 

+j  [+d/c(+u  cos  v+w  sin  v  cos  ,A)~dv  •  w  sin  A+dA  •  u^J 
+k£d#c(v  sin  v  cos  A-u  sin  v  sin  A)+dv(v  sin  A+u  cos  a7] 

and  for  (25) 

£dff’7’]“  +i£-d*(w  sin  05-v  cos  ©  cos  a+da  •  w-do  •  v  sin  aj 

+j  £-d<(u  cos  as  cos  a-w  cos  05  sin  a)+do5(;u  sin  a+w  cos  a)J  « 

+k  j--d/r(v  cos  ©  sin  a-u  Sin  05)-da  *  u-da  .  v  cos  a] 

Consequently,  we  obtain  with  formula  (16)  or  (18)  from  formulas 
(22)  and  (23). 

("r"  +  dr)=  +i£n-dA(w  sin  v  sin  A+v  cos  v)-dv  .  w  cos  A^-dA.  •  v-dd  cos  A  sin  v 
-dx  sin  A-dy  cos  A  cos  y-dx(sih  A  cosytf+cos  A  cos  v  sin*c) 

.  +dy(siri ‘A  sin  a: -cos  A  cos  y  cos^O+dxJ  (26) 

+jJV+d/<,(u  cos  v+w  sin  v  .cos  A)-dy  •  w  •  sin  A+dA  •  u-dd  sin  A.  sin  v, 
+dx  cop  A-dy  sin  A  cos  v+dx(cos  A  cos K  -sin  A  cos  v  sin*-) 
rdy(cps  A  sin/f+sin  A  cos  v  cos>c)+dyL^ 

+kjw-d<r(v  sin  v  cos  A-u  sin  v  sin  A+dv(v  sin  A+u  cos  A)+dd  cos  v  ‘ 

. . -dy  sin  v-dx  sin  v  sin/C-dy  sin  v  cos/c+dzj 

*.4.  *  • 


'  -piJr 
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+3r^  =  +iju-d»c(w  sin  co-v  cos  co  cos  a)+da  ®  w-cko  ®  v  •  sin  a+&d...sin  a  cos  oo 
+dx  cos  a-dy  sin  a  sin  co+dx(cos  a  cos/£+sin  a  sin  tC  sin  co) 

+dy(cos  a  sin>f“sin  a  cosXsin  co-)+dxtJ 

j>v-d^(u:  cos  co  cos  a-w  cos  co  sin  a)+dco(u  sin  a+w  cos  a)+dd  sin  co 

■  -+dy  cos  co-dx  C03  co  sinx+c^  cos  co  cos**+dyJ 

+k|V“=,dX(v.iCos  co  sin  a-u  sin  co)-dct  °  n-dto  .»  v  cos  a+dd  cos  a  cos  co  (27) 
-dx  si.b  a-dy  cos  a  sin  co+dx(-sin  a  cos /C+  cos  a  sin  co  sin/C) 

”dy(sih  a  sin^C+cos  a  sin  co  cos/c  )+dz^ 

If  we. -intersect  the  vector  (  r  +  dr  )  with  a  plane,  we  have: 

X  (r  +  dr)  =  iX(u  +  8y.)+jX(v  +  Sv)+kX(w  +  Sw) 

For  a  plane  at  an  elevation  w  -  constant,  we  obtain, 

■  t-  .  *  < 

X(w  +  $w)  ■  w 


and  with  sufficient  approximation 


w  +.  —  8  w 

v 


X  =  (1  •  i  6  w) 

and,  therefore,  neglecting  second  order  terms: 

X(r”  +  dr)  =  i(u  -  —  6.  w  +  8  u)  +  j(v  -  —  6  w  +  8  v)  +kw  (2C 

w  .  w 

The  point  in  which  the  vector  ~r  +  dr  intersects  the  plane  w  =  constant 
may  be  denoted  by  r'j*’’  and  the  components  of  its  differential  change  . 
within  this  plane  are-  du  and  dve 

Then 


XCr" +  dr)  =  i(u  +du)  +j(v  +  dv)  +  kw 
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and  consequently  with  formula  (28) 


du  *  -  -  J  w  +  8  u 

w 

dv  =  -  -  8  w  +  8  v 

w 


Further,  the  object  vector  R  may  be  expressed  by 
R  a  p.  «  r 

where  p.  is  a  parameter* 

From  Figure  1  we  obtain: 

-*>  A  , 
r  •  l  -  d. 


with  R  •  i  °  |i  •  ? «  5 


we  have 


_  Ri 
■»t  ~  T 


and  with  formulas  (1)  and  (2) 

=  -  X  cos  A.  sin  v  -  Y  sin  A  sin  v  +  Z  cos  v 


-  D 


From  formula  (30),  it  follows  directly  that 

I  =  (1  u 
Y  =  {i  v 
Z  «  p.  w 

Substituting  formulas  (33)  into  formula  (29),  we  obtain 
du  =>  -  ^  ftw  +  %u 

dv  =  -  ^  5w  +  8v 


(29) 


(30) 


(31) 


(32) 


(33) 


(3U) 
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2 

-d  k.  [+  ™  sin  v  sin  A  fr”)  sin  v  cos  A  -  I  cos  v] 


-?dv  |+  ^  co3  A  +  +  zj  sin  aJ 


in  v)] 


♦«'[x]  . 

-dd  £|jl^  003  v  >+  sin  A  sin 
+dx  £  -(jl  cos  aJ 

+dy  ^  sin  v  -  sin  A  cos 

+dx  sin  v  sin  K+  cos  A  cc 

+dy  simy  cos  K  -  cos  A  sir 

_  _  “^diere  «.  x  cos  A  six 

Ml]  ■“ - 


sin  v  sin  K  +  cos  A  cos  K  ~  sin  A  cos 


Y 

n  sirtv  cosK-  cos  A  sinK~  sin  A  cos 


is  v  sin/d)^] 

V  cos^J 


-fthere 


_  -  X  cos  A.  sin  v  -  Y 'Sin  A  sin  v  +  Z  cos  v 


and  the  corresponding  expressions  in  the  tilt  system  are: 

dX  =  +dtf  £“(^-  +  3in  ®  +  ^  cos  co  sin  a  +  Y  cos  co  cos 

tdalr*2J 

+dco  £  ^r-  cos  a  -  T  sin  cQ 

+dd  [  (-  ^  cos  a  cos  co  +  sin  a  cos  co^"J 

+dx  ^  M-  (  ^  sin  “  +  cos  a)J 

+dy  t  p  (  ^  cos  a  sin  co  -  sin  a  sin  co^J 
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+dx.  ^  £sin  a  cos  df -cos  a  sin  co  sink^+cos  a  cos  *+  sin  a  sin  rffsin  <*)] 

(35)* 

+dy  ^  {sin  a  sin>C+cos  a  sin  co  cos/^+cos  a  sintf-  sin  a  cos  K  sin  ca)  J 

Ml]' 


and  ,r 
dT  =■+&*: 


y2  \  YX 

— -  +  Z  1  cos  to  sin  a  sin  co  -  X  cos  co  cos  a1 


[f  *  •) 

pcxV 

U~J 


1  +dx 


EM  cos  a  +.X  sin  a] 
t(-  ^  C03  a  cos  co  +  sin  co^J 
[w-  |  sin  a] 

LU  cos  a  sin  co  +  cos  co^J 


,+<fr 


cos  a  sin  co  +  cos  co 


Of(i{- 

[I] 


sin  a  cos  K-  cos  a  sin 


in  co  sin  -  cos  co  sin>c)J  (39)* 


+dyL  '' 


sin  a  sin  K+  cos  a  siil  co  cos  K  f- +  cos  co  cosXj 


hehe: 

H  =  +  X  3in  a  cos  co  +  Y  sin  co  +  Z  cos  ■  a  co3  co 

d 


*  The  .poefficients  in  formulas  (38)  and  (39)  for  the  .differentials  of  the' 
elements  of  exterior  orientation  agree  with  the  results  obtained  by  4* 
Brandenberger  in  Fehlertheorie  der  iusseren  Orientierung  von  Steilaufnato,- 
meny  Zurich,  191*6..  ■  '  -x  *  s  ; 


With  2  =  1  (unit  distance)  the  formulas  (36),  (37),  (38).,  and  (3 9) 
reduce  to  the  formulas  for  the  differential  changes  d  £  and  d  yj 
of  the  corresponding  standard  coordinates  £  and  »7.* 
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III.  THE  TRIANGULATION  OF  THE  HOST  PROBABLE  POINT  AND  THE  PROPAGATION 
OF  RESIDUAL  ERRORS  IN  THE  TRIANGULATION  PROCEDURE 


In  intersection  photogrammetry,  we  have  in  general  as  many  "rays" 
to  the  target  point  as  there  are  established  measuring. stations.  Due 
to  the  inevitable  errors,  these  rays  will  not  intersect  in  space  unless 
corrections  are  applied  to  the  individual  rays.  The  criterion  selected 
for  the  most  probable  corrections  must  lead  to  the  determination  of  the 
most  probable  point  of  intersection.  First  we  shall  assume  that  the 
systematic  errors  of  the  orientation  elements  are  small  compared  with 
the  expected  accuracy  of  the  plate  measurements  and  can,  therefore,  be 
neglected.  This  means  that  We  assume  that  the  elements  of  interior  and 
exterior  orientation  are  known  sufficiently  a  ccurately  either  from  dial 
rea'dings  and  instrument  calibration  procedures,  or  from  preceding  compu¬ 
tations  of  the  orientation  elements  from  recorded  control  points.  Conse¬ 
quently,  we  may  say  that  the  most  probable  point  is  determined  by  the 
intersection  of  such  rays  for  which  the  corrections  applied  to  the  measured 
-  plate  coordinates  of  the  individual  measuring  stations  satisfy  the  condi¬ 
tion  that  the  sum  of  the  least  squares  is  a  minimum.  After  the  most 
probable  corrections  have  been  applied  we  have  a  rigorous  triangulation 
and  the  spatial  coordinates  may  now  be  computed  from  any  two  stations. 

Formulas  for  such  a  triangulation  follows 


In  Figure  5  we  have  two  stations.  -  I  and1  II  -  for  which  -  b  +  P^ » 
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With 


P1  =  iul  +  ^V1  +  lcwl 
=  iu2  +  jv2  +  kw2, 


b  =  ib  +  ib  +  kb_ 
x  J  y  2 


the  vector  products  aret 

i 


b  x  p2  - 


0 

b 


u„ 


P^_  ?  P2 


u. 


u„ 


k 

b 

2 

w. 


W-, 


*  2  **  2 

The  directions  of  both  vector  products  are  the  same  and  consequently 


|J-  = 


L 

b 

b 

b. 

b 

b 

b 

X 

7 

7 

z 

X 

-  z 

u. 

V. 

v„ 

U- 

w„ 

2 

2 

2 

.  2 

2 

2 

*1 

V1 

V1  • 

W1 

*1 

W1 

U2 

V2 

V2. 

CVJ 

J2 

U2 

W2 

Oil) 


Directly  from  formula  (LtO)  : 
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x  =  \  "  P2  U2 


Y  =  V1  =  ^2  V2  +  b 
Z  “  ^  W1  “  ^2  W2  +  \ 


The  u,  v  and  w  values  are  conputed  from  formulas  (3)  or  (9).  It  is 
noted  that  the  vector  triplet  of  the  Station  II  must  first  be  trans¬ 
formed  in  the  vector  triplet  of.  the  Station  I  in  order  to  eliminate 
the  earth's  curvature „ 

Although  the  result  obtained  with  (^2)  is  correct  only' for  a 
unique  solution,  these  formulas  enable  us  to  study  the  propagation 
of  the  error  of  the  individual  stations  in  a  triangulation  proce¬ 
dure,,  We  shall  limit  our  consideration  to  two  stations.  Further, 
we  shall  neglect  the  influence  of  the  curvature  of  the  earth  and 
consider  the  errors  on  the  geodetic  basic  data  as  negligibly  small. 

We  introduce  a  cartesian  coordinate  system  with  its  drigin  at  Station  ,1 
oriented  in  such  a  way  that  the  horizontal  base  line  component  coin¬ 
cides  with  the  +  y  axis,  as  shown  in  Figure  6. 


** - ° - 1  ^,6 

From  formulas  Xli?),  we  obtain  by  differentiation,  the  differential  changes 
in  the  space  coordinates,  8X,  &Y  and  8Z. 

&X  =  p,  8  u-  +  u,  8  LU 


&Y  =  p1  8  vL  +  v^  S  p-^ 

Pi  $  «i  +  Wi  &  Pi  +  p2  8  w2  +  w2  8  p2 
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With  formulas  (Ul)5  (h.2)  and  (1*3) 

T|alSul"  ^Y"b  ^  t12&u2~X  ^  |1l8vl"(J'2Sv2  ^ 

M - 5 

I^r£i(^1»^26'>2)-(Y-b)^1''^26v2 

81  b  (UU) 

Z(2T~b)-bzY  2Z-bz 

6Z  _  2bX  ^ ”d ^  ""  2b~^^l8vl"^i2^V2  )+,2  ^i^wi+iJ,28w2 ^ 


where  Sin,  8v  and  &w^  are  obtained  from  formulas  (26)  or  (2li) 
From  formulas  (i±2)  and  (26) P  we  have; 

u.&u.  =  -Y.  «dA. 
t.  i  i  i . 


-Z.  cos  A.dV. 
i  ii 


"(Z^  sin  sin  A^+Y^  cos  v^)»d4^ 

-ji^  sin  A^«dx^ 

-u.  cos  A.  cos  v.«dy. 

T.  i  i  J  i 

-{i^  (sin  cos  /c^+cos  A^  cos  vi  sin  *\  )d!x^ 

t 

+ll.  (sin  A.  sin  K.  -cos  A.  cos  v.  cos  tr.  s 

^l  i  1  l  l  ^i)dy. 


~ll;  cos  A,  sin  v.  »dd. 

-  l  i  l  l 


ii-6v.  =  +XdA. 

l  i  l 


-Z  sin  A. •dv. 

•ii 

+  (X  cos  v.+Z.  sin  v.  cos  A.  )*d/<. 

ii  i  i  i 

+  ll.  cos  A'.  »dx. 

^i  ii 

~|i.  sin  A.  cos  v.  dy. 

X  1  X  1 
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+LL-.  (cos  A.  cos  k  ■  "Sin  A.  cos  v.  sin*.)'"dx- 
i.  i  i  i  1  i  x. 

-u.fcos  A.  sin*.+sin  A.  cos  v.  cos  K , )  «dy. 

^l'  x  x  i.  i  i'  J  x 

-u.  sin  A.  sin  v.  ®dd. 
ri  x  ix 

*J'i&wi-'S  +  (Y^  sin  A,  t  X  cos  A-^)  •dvi 

/ 

-(Y.  cos  A.  sin  v.~X.  sin  A.  sin  v.)«d*.  . 

x  x  i  i  i  l  l 

sin  vi«dyi 


*^i 

sin 

vi 

sin  K-  8dx. 
i  -  i 

■^i. 

sin 

vi 

cos 

cos 

vi 

«dd. 

i 

where,  in  accordance  with  formulas  (31)  and  (32), 

Di 

a.  “  -r-  and  D.  ~  -X  cos  A.  sin  v.  -I.  sin  A.  sin  v. 

'i  d.  x  ■  ,  1  x  xi  l  i 

and  from  formula  (27) 

-Y^  sin  cu  •dca^ 

-(Z.  sin  a. “Y.  cos  co.  cos  . a.  )*d*. 
v  x  xi  x  i'  i 

cos  a^»dx^ 

j 

sin  sin  ecu  •  dy\ 


/ 


+ 


+a.  (cos  a.  cos/c4+sin  a.  sinK’. 

^i  '  i  i  i  i. 

+ja^  (cos  sinAC^-sin  cosAc^ 


sin  cd^)  •  dx^ 
sin  gk  )  •  djj!\ 


cos 
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+i_i^  sin  cos  ox  •  dcL 

u.&v.  =  +  (X  sin  a.  +  2.  cos  a.  )dco. 

•l'  1  •  • '  •  •  i  i  i,i 

.  * 

-(X  cos  co^  cos  -  Z^  cos  co^  sin  a^)«  d/c^ 

'  +ii.  cos  co.  •  dy. 

i  i  i  • 

-a-  cos  co.  sin/c."  dx. 

t.  ill 

+ii^  cos  co^  cos/c^»  dy^ 

+ii.  sin  co.  ®  dd. 

T.  i  i 

P-i8wi  °  "Xcet^ 

-T^  cos  dco^ 

-(Y.  cos  co.  sin  a.  -  X.  sin  co. )®  d/C. 

i  i  i  i  i  i 

“(i^  sin  a^o  dbci 

/ 

-ll.  cos  a.  sin  co.  dy. 

^i  l, ,  .  i  ^i 

-a. (sin  a.  cos /C. -cos  a.  sin.  co.  sin  K- )  •  c£x. 

i  x  i  x 

-a.  (sin  a.  sin  >f.+co’s  a.  sin  co.  cos>e.)»  dy. 

pi  i  •  i  i.-i.  i,  l 

*  \ 

cos  cos  co^«  dd^  . 

where,  with  formulas.  (31)  and  (32), 

Di 

a.  c  3—  and  D.  =  X  sin  .a,  cos  co.  +  Y.  sin  co.  +  Z.  cos  a.  cos'  co. 
yi  d.  i  .  i  l.  1  l  i.  i  i . 

Oiu^'CL, 

Helljirut  Schmid 
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